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ABSTRACT 

Aims. We present a spectral analysis of a deep (220 ks) XMM-Newton observation of the Phoenix cluster (SPT-CL J2344-4243). We 
also use Chandra archival ACIS-I data that are useful for modeling the properties of the central bright active galactic nucleus and 
global intracluster medium. 

Methods. We extracted CCD and reflection grating spectrometer (RCS) X-ray spectra from the core region to search for the signature 
of cold gas and to finally constrain the mass deposition rate in the cooling flow that is thought to be responsible for the massive star 
formation episode observed in the brightest cluster galaxy (BCC). 

Results. We find an average mass-deposition rate of M = 620 (-190 -l- 200),,,(-50 150)jj.j,Mo yr“* in the temperature range 

0.3-3.0 keV from MOS data. A temperature-resolved analysis shows that a significant amount of gas is deposited at about 1.8 keV 
and above, while only upper limits on the order of hundreds of Mq yr“' can be placed in the 0.3-1.8 keV temperature range. From pn 
data we obtain M = 210 (-80 + 85)j,o, (-35 -l- 60 )j3,j,Mq yr“’ in the 0.3 - 3.0 keV temperature range, while the upper limits from 
the temperature-resolved analysis are typically a factor of 3 lower than MOS data. No line emission from ionization states below Fe 
XXIII is seen above 12A in the RCS spectrum, and the amount of gas cooling below ~ 3 keV has a formal best-fit value M = 122/|22 
Mq yr“'. In addition, our analysis of the far-infrared spectral energy distribution of the BCC based on Herschel data provides a star 
formation rate (SFR) equal to 530 Mq yr“' with an uncertainty of 10%, which is lower than previous estimates by a factor 1.5. Overall, 
current limits on the mass deposition rate from MOS data are consistent with the SFR observed in the BCG, while pn data prefer a 
lower value of M ~ S FRI3, which is inconsistent with the SFR at the 3cr confidence level. 

Conclusions. Current data are able to firmly identify a substantial amount of cooling gas only above 1.8 keV in the core of the Phoenix 
cluster. At lower temperatures, the upper limits on M from MOS and pn data differ by a factor of 3. While the MOS data analysis 
is consistent with values as high && M ~ 1000 within 1 cr, pn data provide M < 5Q0Mq yr“* at 3a confidence level at a temperature 
below 1.8 keV. At present, this discrepancy cannot be explained on the basis of known calibration uncertainties or other sources of 
statistical noise. 
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1. Introduction 

The majority of baryo ns in clusters of galaxies is constituted 
by virialized hot gas (iLin et al.l I2003I : iGonzalez et al.l l2013l) 
that emits X-ray via thermal bremsstrahlung. Temperature, den¬ 
sity, and chemical composition of the so-called intracluster 
medium (ICM) can be directly measured thanks to X-ray imag¬ 
ing and spectroscopic observations. Spatially resolved spectro¬ 
scopic studies showed a significant temperature decrease and 
strongly peaked surface brightness profiles in the center of a 
significant fraction of the cluster population. The short cool¬ 
ing times associated with these high-density gas regions led to 
the conclusion that a massive cooling flow was developing in 


the ICM in most of the clus ters (| Silk j_97^ Cowe^^BinnevI 
[T97llK^ & Nulsenl[l97^lMathews & Bregmanlll978l) . The 


fate of this cooling gas would be to feed massive star formation 
episodes. 


On the bas is of the isobaric cooling-flow model 
(iFabian & NulsenI I1977I : iFabianI Il994ll . it was estimated 
that typical cooling flows may develop mass deposition rates 
in the range of a few xlOO - IOOOMq yr *. However, the lack 
of massive star formation events and of large reservoirs of cold 
gas at the center of galaxy clusters casts some doubts on the 
hypothesis of a complete cooling of the ICM in cluster cores. 
The picture changed dramatically when X-ray observations, in 
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particular of the RGS instrument onboard the XMM-Newton 
telescope, revealed a severe deficit of emission lines compared 
to the predictions of the isobaric cooling-flow model in all 
the groups and clusters of galaxies with putative cooling flows 


jTamurae^^ 


iPeterson et al 


12001: IPeterson et al.l l200lt iKaastra et alJ l200ll 


2003h . Interestingly, this result has also been 


found in XMM-Newton and Chandra CCD spectra, despite 
the lower resol ution, thanks to the prominent complex of iron 
L-shell lines ( McNamara et alJ 2 000t Bohringeret^ 1200 It 


iMolendi & Pizzolatoll200lT Bohringer et al.ll2002h . This directly 
implies that the cooling gas is present only in small amounts, 
typically ten times lowe r than expected for steady- state, isobaric 
radiative cooling (see IPeterson & Fabiaiil l2006l) . This deter¬ 
mined a change from the cooling-flow paradigm, with typical 
deposition rates of about 100 - IOOOMq yr"\ to the cool-core 
paradigm, where most of the gas is kept at temperatures higher 
than one-third of the ambient cluster temperature, and the mass 
deposition rate, if any, is due to a residual cooling flow of about 
few tens Mq yr"'. 

Another direct implication of these observations is that there 
must be some process that heats the gas and prevents its cool¬ 
ing. Among the many mechanisms investigated in the past 
years, feedback from active galactic nuclei (AGN) is consid¬ 
ered the most plausible heating sou rce. Radio AGN are ubiq¬ 
uitous in cool cores (see ISunll2009h . and interactions between 
the radio jets and the ICM have been observed unambigu¬ 
ously. AGN outb ursts can in principle in ject sufficient energy 
into the ICM (see lMcNamara et'^l2005h . The relativistic elec¬ 
trons in jets associated with the central cluster galaxy are able 
to carve large cavities into the ICM. The free energy associ¬ 
ated with these bubbles is plausibly tra nsferred into the ICM 
and thermalized through turb ulence (see iMcNamara & NulsenI 
1201 2t IZhuravleva et alJF2014l) . In addition, there is increasing 
evidence of interactions of AGN outflows with metal-rich gas 
along the cavities and edges of radio jets of some individua l 
clusters and groups (iKirkpatrick et al.ll201 iHEttori et al.ll2013bll . 
consistent with numerical simulations showing that AGN out¬ 
flows are able to advect ambient, ir on-rich material from the core 
to a few hundred kpc away (e.g.. iGaspari et al]l201 lall5) . The 
feedback mechanism has been observed in its full complexity i n 
nearby clusters such as Perseus ( Fabian et al.ll2003ll20()6ll201 Ih . 
Hydra A (iMcNa mara et alj|2000ll . and a few other clusters (see 
iBlanton et al.ll20lIlK 

In addition, the regular behavior of cool cores is not only 
observed in local clusters, but seems to hold up to high red- 
shifts. High angular resoluti on observations of cool cores at 
z < 1 performed by our group (ISantos et al.ll201(il2012h showed 
that the radio feedback mechanism is already present at z ~ 1. 
Temperature and metallicity profiles in cool cores are broadly 
consistent with local ones, with a remarkable difference in the 
the metal distribution that ap pears to be more conc entrated in 
the core than in local clusters (IDe Grand! et al.ll20l4) . This indi¬ 
cates that radio feedback also plays a role in the spatial distribu¬ 
tion of metals. The regularity of the cool-core appearance over 
the entire cluster population and a wide range of epochs points 
toward a gentle heating mechanism, with the radio AGN acting 
with a short duty-cycle to counterbalance the onset of cooling 
flows since the very first stages of cluster formation. At the same 
time, outburst shocks may provide a more violent heating mech¬ 
anism. It is poorly understood whether AGN heating of the ICM 
occurs violently through shocks or through bubbles in pressure 
equilibrium that cause turbulence. Outburst shock with jumps in 
temperature are very hard to detect, and shocks have been un¬ 
ambiguously detected in only a few cases (see, e.g., A2052 and 


NGC.5813. iBlanton et al.llMH iRandall et al.l 1201 ll) . It is now 
widely accepted, however, that the mechanical energy provided 
by the AGN through jets is sufficient to overcome the cooling 
process in cluster cores. Therefore questions remain about the 
detailed physical mechanism that transfers energy to the ICM, 
and how this mechanism gives rise to the regular cool-core ther¬ 
mal structure, which requires a minimum temperature a factor 
~ 3 lower than the ambient cluster temperature. 

Surprisingly, a recent observation introduced further 
changes in the picture outlined here. The SZ-selected clus- 
ter SPT-CLJ2344-424 3 (also known as the Phoenix cluster, 
iMcDonald et al.ll2012l) at z ~ 0.596 for the first time shows 
hints of a massive cooling-flow-induced starburst, suggesting 
that the feedback source responsible for preventing runaway 
cooling may not yet be fully established. SPT-CLJ2344-4243 
shows a strong cool core with a potential mass deposition 
rate of ~ 3000Mpi yr ~* derived from the X-ray luminosity. 
iMcDonald et al.l (1201 2l) argued that the Phoenix cluster might 
harbor an almost isobaric cooling flow with an unusually high 
mass deposition rate. The strongest hint comes from the very 
high star formation rate (SFR) observed in the brightest clus¬ 
ter galaxy (BCG), which has originally been estimated to be 
~ 70 0Mo vr~*, with large Icr errors ranging from 200 to 5QQMq 
yr*' (IMcDonald et al.l [2012h . However, an accurate measure¬ 
ment of the SFR is made difficult by the presence of a strongly 
absorbed AGN, whose contribution to the BCG emission can be 
accounted for in differe nt ways. Recently, the H STAVFC3 obser¬ 
vation of the Phoenix (IMcDonald et al.ll2013h showed filamen¬ 
tary blue emission out to 40 kpc and beyond, and the estimated, 
extinction-corrected SFR has been updated to a more accurate 
value of 798 + 42Mo yr“', consistent with optical and IR data at 
lower spatial resolution. 

In this work we present the analysis of a 220 ks observation 
with XMM-Newton awarded in AO 12 on the Phoenix cluster, 
with the main goal of investigating the thermal structure of the 
cool core and comparing the mass deposition rate in the core to 
the star formation rate in the BCG. We also use archival Chandra 
data (about 10 ks with ACIS-I) to model the emission of the 
AGN in the BCG and global ICM properties. Finally, we revise 
the SFR in the BCG on the basis of far-IR (FIR) data from the 
Herschel Observatory. 

The paper is organized as follows. In Sect.|2]we describe the 
reduction of the XMM-Newton and Chandra data. In Sect. Owe 
present the results from the Chandra data analysis on the central 
AGN spectrum and the global ICM properties. In Sect. |4]we de¬ 
scribe our analysis strategy and present the results on the cool- 
core temperature structure from EPIC MOS and pn data, RGS 
data, and Chandra data. In Sect. O we revise the measurement 
of the star formation rate in the BCG from FIR data in view of 
a comparison with the mass deposition rate. Finally, our conclu¬ 
sions are summarized in Sect.|6] Throughout the paper, we adopt 
the seven-year WMAP cosmolog y with Da - 0-73, = 0.27, 

and Ho - 70.4 km s ' Mpc * dkomatsu et al.ll20lTl) . Quoted 
errors and upper limits always correspond to a 1 cr confidence 
level, unless stated otherwise. 

2. Data reduction 

2.1. XMM-Newton: EPIC data 

We obtained a total of 225 ks with XMM-Newton on the 
Phoenix cluster in A012Q. Data were acquired in November and 

' Proposal ID 72270, ’’The thermal structure of the cool core in the 
Phoenix cluster”, PI P. Tozzi 
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Table 1. XMM-Newton data: exposure times for each Obsid after data reduction. 


Obsid 

EPIC detector 

effective texp 
ks 

0722700101 

MOSl 

128.0 

0722700101 

MOS2 

128.0 

0722700101 

pn 

103.0 

0722700201 

MOSl 

92.0 

0722700201 

MOS2 

92.0 

0722700201 

pn 

81.5 

0693661801 

MOSl 

16.0 

0693661801 

MOS2 

16.5 

0693661801 

pn 

8.6 


December 2013 (Obsid 0722700101, 132 ks, and 0722700201, 
93 ks). We added a shorter 20 ks archival observation taken in 
2012 to our analysis. 0 

The observation data files (ODF) were processed to pro¬ 
duce calibrated event files using the most recent release of the 
XMM-Newton Science Analysis System (SAS vl4.0.0), with 
the calibration release XMM-CCF-REL-323, and running the 
tasks EPPROC and EMPROC for the pn and MOS, respectively, 
to generate calibrated and concatenated EPIC event lists. Then, 
we filtered EPIC event lists for bad pixels, bad columns, cosmic- 
ray events outside the field of view (EOV), photons in the gaps 
(PLAG=0), and applied standard grade selection, corresponding 
to PATTERN < 12 for MOS and PATTERN <= 4 for pn. We 
removed soft proton flares by applying a threshold on the count 
rate in the 10-12 keV energy band. To define low-background 
intervals, we used the condition RATE < 0.35 for MOS and 
RATE < 0.4 for pn. We found that for Obsid 0722700101, we 
have about 128 ks for MOSl and MOS2 and 103 ks for pn after 
data reduction. Eor Obsid 0722700201 we have 92 ks for MOSl 
and MOS2 and 81.5 ks for pn. This means that removing high- 
background intervals reduces the effective total time from 225 ks 
to 220 for MOS (a loss of 2% of the total time) and to 184.5 ks 
for pn (a loss of 18% of the total time). The archival data, Obsid 
0693661801, were affected by flares by a larger amount. Erom 
20 ks of exposure in Obsid 0693661801, we have 16 and 16.5 
ks for MOSl and MOS2, respectively (a loss of about 18-20%), 
and 8.6 ks for pn (a loss of about 60% of the total time). 

Eor each Obsid we merged the event files MOS 1 and MOS2 
to create a single event file. This procedure has been adopted 
because the MOS cameras typically yield mutually consistent 
fluxes over their whole energy bandpass 0. Einally, we also re¬ 
moved out-of-time events from the pn event file and spectra. In 
Table 1, we list the resulting clean exposure times for the pn and 
the MOS detectors for the three exposures used in this work. 

Effective area and response matrix files were generated with 
the tasks arfgen and rmfgen, respectively, for each obsid. Eor 
MOS, effective area and response matrix files were computed for 
each detector and were eventually summed with a weight corre¬ 
sponding to the effective exposure time of MOSl and MOS2. 


^ Proposal ID 069366, PI M. Arnaud. 

^ See http://xmm2.esac.esa.int/docs/documents/CAL-TN-00I8.pdf 


2.2. XMM-Newton: RGS data 

We reduced the RGS data set using the standard SAS vl4.0.0 
pipeline processing through the RGSPROC tool. We filtered soft 
proton flares by excluding time periods where the count rate 
on CCDS|3 is lower than 0.1 cts^ in a region free of source 
emission. The resulting effective exposure times, after remov¬ 
ing flares from the data, are listed in Table |2 As background we 
adopted the model background spectrum created by the SAS task 
RGSBKGMODEL, which can be applied to a given observation 
from a combination of observations of empty fields, based on the 
count rate of the off-axis source-free region of CCD9. We veri¬ 
fied that the background spectrum obtained in this way is entirely 
consistent with a local background extracted from beyond 98% 
of the RGS point spread function (PSE). Einally, we focused on 
the first-order spectra and combined spectra, backgrounds, and 
responses from all observations and from both RGS instruments 
using the SAS task RGSCOMBINE. 

Table 2. Cleaned exposure times for RGS detectors 


Obsid 

RGS detector 

effective texp 
ks 

0722700101 

RGSl 

126.4 

0722700101 

RGS2 

125.2 

0722700201 

RGSl 

92.4 

0722700201 

RGS2 

92.7 


2.3. Archival Chandra data 

The Phoenix cluster has been observed with ACIS-I for 11.9 ks 
in the VEAINT mode (Obsid 1340lfl We performed a standard 
data reduction starting from the levels 1 event files, using the 
CIAO 4.6 software package, with the most recent version of the 

This is the CCD that generally records the fewest source events be¬ 
cause of its location close to the optical axis and is the most susceptible 
to proton events. 

^ Proposal ID 13800933, ’’Chandra Observation of the Most Massive 
Galaxy Clusters Detected in the South Pole Telescope Survey”, PI. G. 
Garmire. 
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Chandra Calibration Database (CALDB 4.6.3). We ran the task 
acis_process_events to flag background events that are most 
likely associated with cosmic rays and removed them. With this 
procedure, the ACIS particle background can be significantly re¬ 
duced compared to the standard grade selection. The data were 
filtered to include only the standard event grades 0, 2, 3, 4, and 
6. We visually checked for hot columns that were left after the 
standard reduction. As expected for exposures taken in VFAINT 
mode, we did not find hot columns or flickering pixels after fil¬ 
tering out bad events. Finally, we filtered time intervals with high 
background by performing a 3cr clipping of the background level 
using the script analyze_ltcrv. Only a negligible fraction of 
the exposure time was lost in this step, and the final effective ex¬ 
posure time is 11.7 ks. We remark that our spectral analysis will 
not be affected by possible undetected flares, since we are able 
to compute the background in the same observation from a large 
source-free region close to the cluster position, thus taking into 
account any possible spectral distortion of the background itself 
induced by the flares. 

3. Spectral analysis of Chandra data: AGN spectrum 
and global ICM properties 

We first performed the spectral analysis of the Chandra data. 
Despite the short exposure tim43, the high angular resolution of 
Chandra can provide important parameters useful in the spectral 
analysis of the XMM-Newton data. In particular, it is immedi¬ 
ately clear from the high-resolution hard-band image that the 
BCG of the Phoenix cluster hosts a powerful obscured AGN. 
In XMM-Newton data, the point spread function has an half¬ 
energy width (HEW) of about 15” at the aimpoint, which causes 
the AGN emission to be spread over the entire starburst region. 
Chandra data allow us to accurately measure the spectrum of the 
central AGN and eventually model its emission in the analysis 
of the XMM data. 

We extracted a circular region with a radius of 1.5 arcsec 
to analyze the position of the central source in the hard-band 
Chandra image (RA = 23:44:43.9, DEC = -42:43:12.64). At 
variance with most cool core clusters, the AGN in the center of 
the Phoenix is extremely X-ray luminous. We only considered 
the energy range 1.0-10 keV for fitting purposes to avoid resid¬ 
ual contamination from the thermal emission of the ICM. We 
detected 575 net counts in the 1.0-7 keV band, most of them in 
the hard 2-7 keV band. In ou r fit we fixed the redsh ift to the 
optical value Zopt - 0.596 (see lMcDonald et^l2012h . We also 
fixed the Galactic absorption to the va lue NHcmi - 1.52 x 1 0^^ 
cm“^ obtained from the radio map of iKalberla et alJ (1200.51) at 
the position of the cluster. The AGN was modeled as a power 
law with an intrinsic absorption (XSPEC model zwabs x pow) 
convolved by the Galactic absorption model (tbabs). The slope 
of the power law was frozen to the value F = 1.8. 

We found and intrinsic absorption NH = (46 + 7) x 10^“ 
cm“^ and unabsorbed intrinsic luminosities of (3.1 + 0.1) x 10"^^ 
and (4.8 + 0.2) x 10"^^ erg s^* in the 0.5-2 keV and 2-10 keV 
bands, respectively. The Chandra spectrum of the AGN and the 
best-fit model are shown in Fig. |2] According to this, the cen¬ 
tral AGN can be classified as a type II QSO. T hese results 
are consistent with the results of lUeda et al.l (l2013l) . who found 
NH = (32 + 9) X 10^^ cm“^ and F = 1.54 + 0.27 from the 
combined analysis of Suzaku XIS and HXD and Chandra (with 
CALDB 4.5.3). We also note that if we left F free in our fit. 


® A new 100 ks observation with Chandra has been taken in August 
2014, PI M. McDonald. 


we found a much flatter spectrum with an intrinsic absorption 
lower by a facto r of two, in agreement with the findings of 
lUeda et 5] (1201 3l) for the Chandra data alone. This is due to the 
well-known degeneracy between spectral slope and intrinsic ab¬ 
sorption. Given the very low value ~ 0.6 for the intrinsic spectral 
slope obtained in this way, we prefer to rely on the results with 
F = 1.8 consistent with the Suzaku+Chandra analysis. Since the 
contribution of the AGN in the soft band is crucial in our spec¬ 
tral analysis of the core region, we will eventually allow NH to 
range from 23 to 53 x 10^^ cm“^, to span the upper and lower 
Icr limits of the two measurements in the analysis of the XMM 
data. Finally, we note that we cannot clearly identify the neutral 
Fe line at 6.4 keV rest-frame when adding an unresolved line 
component. We also remark that our analysis of Chandra data 
was obtained by removing the surrounding ICM emission, and 
n ot by model i ng the thermal and AGN components together as 
in lUeda et al.l (120131) . 

The total ICM emission contributes about 6400 net counts in 
the 0.5-7 keV within a radius of about 650 kpc, beyond which 
the surface brightness reaches the background level. This en¬ 
abled us to perform the spectral analysis in independent rings 
with slightly fewer than 1000 net counts each. The projected 
temperature and iron abundance profiles are shown in Fig. [T]in 
the left and right panels, respectively. The temperature profile 
clearly shows a prominent cool core with a decrease of at least 
a factor of 2 from 200 kpc to the inner 50 kpc. This is mirrored 
in the iron abundance profile as a clear peak toward the center, 
where the iron abundance rea ches solar metallicitv (with solar 
metallicity values measured bv lAsnlund et al.ll2005h with an un¬ 
certainty of 30%. These results were used to complement the 
spectral analysis of XMM-Newton data. 

We also computed the deprojecte d temperature and densit y 
profiles with a backward method (see lEttori et al.ll^ 1 Oi l20 1 3^ . 
which makes use of the geometrically deprojected X-ray sur¬ 
face brightness and temperature profiles to reconstruct the hy¬ 
drostatic mass profile. This is assumed to be described by the 
NEW functional form (iNavarro et alJIlOO"^ . The best-fit param¬ 
eters of the mass profile were obtained through the minimiza¬ 
tion of the statistics defined as the sum of the squared dif¬ 
ferences between the observed temperature and the temperature 
estimates obtained by inverting the equation of the hydrostatic 
equilibrium, weighted by the observational errors on the spec¬ 
troscopic temperature. The gas density was derived from the 
deprojected surface brightness and the total mass model. With 
this method we measured a total mass, under the assumption 
of hydrostatic equilibrium, of M^qq — (2.34 + 0.71) x 10'^ Mq 
at Rsoq - 1627 + 235 kpc. The ICM mass within R^oo is 
Micm - (2.1 ± 0.2) X 10*^ Mq. Erom this, we obtained an ICM 
fraction of //cm = 0.09 + 0.03 at R 50 Q. As a simple check, we 
estimated the parameter Tx = (27.3+4.1)xlO*^keV M© assum¬ 
ing as representative of the entire cluster the temperature mea¬ 
sured between 100 and 500 kpc, and compu ted the mass e s timate 
from Yx assuming the relation described in lArnaud et al.l (1201 Ol) 
with the slope fixed to the self-similar value. We found a mass 
of (2.15 + 0.24) X IO'^Mq, consistent with the value previously 
computed from the hydrostatic equilibrium well within 1 cr. 

If we extrapolate the measured profile up to the vir ial radius 
adopting the best-fit NEW profile dNavarro et al.lll99^ . we find 
M 200 = 3.4 + 1.2 X 10*^ Mq at /?200 = 2500 + 400 kpc. We note 
that this value is about a factor of 2 higher than the mass estimate 
M 200 SZ - (1.66 + 0.23stat ± 0.44,vtf) x 10*^ Mq obt ained from 
SZ-mass scaling relations in IWilliamson et akl (1201 ll) . However, 
considering that SZ-inferred masses are found to be statistically 
lower by a factor of 0.78 in their sample, the discrepancy be- 
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Fig. 1. Left panel: projected temperature profile of SPT-CLJ2344 obtained from the 11.7 ks Chandra ACIS-I observation, after 
removing the central AGN. Errorbars correspond to 1 cr. Right panel: projected iron abundance profile of SPT-CLJ2344 from 
Chandra data. 



Energy (keV) 

Fig. 2. Folded spectrum with best-fit model for the AGN in the center of the BCG of the Phoenix cluster. The spectrum is extracted 
from a radius of only 1.5 arcsec centered on the peak of the hard X-ray emission. 


tween X-ray and SZ mass for the Phoenix is reduced to ~ Icr, 
and therefore is not significant. Finally, by extrapolating the ICM 
density distribution, we find M/cm = (3.2 + 0.3) x 10'"^ M© for a 
virial ICM fraction of //cm = 0.095 + 0.035. 


4. Spectral analysis of XMM-Newton data 

4.1. Analysis strategy 

We extracted the spectra from a circular region with a radius of 
13.5 arcsec. This region was chosen to include the bulk of the 
emission from the core region, which is estimated to be con¬ 
fined within a radius of 40 kpc (about 6 arcsec at z ~ 0.6), 
plus 7.5 arcsec corresponding to the half-power diameter of the 
XMM-Newton PSF Since the cold gas is expected to be concen¬ 


trated toward the center of the cluster, we assumed that the signal 
from the cold gas emission outside this radius is negligible. This 
assumption is consistent with the results we obtained from the 
RGS spectra in two different extraction regions (see Sect. 4.3). 

The background was sampled from a nearby region on the 
same CCD that was free from the cluster emission and was sub¬ 
tracted from the source spectrum. We note that the total back¬ 
ground expected in the source region, computed by geometri¬ 
cally rescaling the sampled background to the source area, only 
amounts to 0.3% of the total observed emission. Before perform¬ 
ing the final spectral analysis, we ran a few spectral fits arti¬ 
ficially enhancing the background by factors of the order 1.2- 
1.5, which range encompasses any possible uncertainty in the 
background level. We found no relevant differences in the best- 
fit values as a function of the background rescaling factor. We 
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conclude that the background only weakly affects the spectral 
analysis. We therefore only consider the results obtained with the 
background sampled from the data and geometrically rescaled to 
the source region. 

In XMM-Newton spectra the emission from the central AGN 
is mixed with the ICM emission, therefore we modeled its con¬ 
tribution with a power-law spectrum with a fixed slope of F = 
1.8 and an intrinsic absorption ranging from 23 to 53 x 10^“ 
cm“^ at the redshift of the cluster, as found in Sect. a We 
also left the normalization free for each separate spectrum to 
account for possible differences in the calibration of Chandra 
and XMM-Newton. We consistently found values in the range 
(1.0-1.4)xl0“^ for the normalization of the power-law compo¬ 
nent. The Galactic absorption, instead, was frozen to the value 
NHr.ni - 1.5 2 X 10 ^*^ cm“^ obtained from the radio map of 
iKalberla et al.l (l2005h at the position of the cluster, and it is de¬ 
scribed by the tbabs model. The effects of possible uncertain¬ 
ties on the Galactic absorption are discussed together with other 
systematic effects in Sect. 4.4. 

We assumed that the gas in the residual cooling flow can 
be described by an isobaric cooling model, therefore we use d 
the mkcflow spectral model (iMushotzkv & Szvmkowial3ll988h . 
This model assumes a unique mass-deposition rate throughout 
the entire temperature range T„ii„-Tmax- Since we are interested 
in the gas that is completely cooling and is continuously dis¬ 
tributed across the entire temperature range, we set I'm,,, = 0.3 
ke'V and T„tax - 3.0 ke'V. We set the lowest temperature to 
0.3 ke'V because this is the lowest value that can possibly con¬ 
tribute to the emission in the Chandra and XMM-Newton en¬ 
ergy range, given the relatively high redshift of the Phoenix clus¬ 
ter. However, a single mkcflow model may not be sufhcient to 
investigate the structure of the cool core. The actual situation 
may be more complex, with some of the gas, above a given 
temperature threshold, cooling at a relatively high rate consis¬ 
tently with the isobaric cooling-flow model, while colder gas 
may have a much lower mass-deposition rate. Grating spectra 
of cool cores are traditionally fitted with an isobaric cooling- 
flow mo del with a cutoff tempera ture below which no gas is 
detected (iPeterson & Fabianll2006h . To explore a more complex 
scenario, we aimed at separately measuring the cooling rate (i.e., 
the mass deposition rate) in several temperature bins. In this 
case, the temperature intervals were fixed to 0.3-0.45, 0.45-0.9, 
0.9-1.8, and 1.8-3.0 ke'V. This choice is similar to using the 
single-temperature mekal model for a discrete set of tempera¬ 
tures, but with the advantage of a continuous description of the 
gas instead of a set of discrete temperature values. 

The contribution of the much hotter ICM component along 
the line of sight was modeled with a mekal model with a free 
temperature. Above 3 keV, a single-temperature mekal model 
can account for several hot components because it is not possible 
to resolve the te mperature structure ab ove 3 keV with the spec¬ 
tral analysis (see lMazzotta et afl2004l) . This means that the pos¬ 
sible contributions to the emission from temperatures between 
3.0 keV (the highest temperature of the mkcflow) and the best- 
fit temperature of the mekal model are already described by 
a single-temperature mekal component. Therefore, despite the 
strong temperature gradient toward the center, the relevant part 
of the thermal structure of the ICM is properly treated by assum¬ 
ing a single temperature for the hot component and exploring 
the low-temperature regime with a multi-temperature mkcflow 
model. 

In practice, our fitting method consists of two measurements 
of the mass deposition rates, using the following models: 


- a single cooling-flow model mkcflow plus one single¬ 
temperature mekal component. The lowest temperature of 
the mkcflow component is frozen to 0.3 ke'V, and the largest 
to 3.0 ke'V. By setting the lowest temperature to 0.3 keV, we 
can interpret the normalization of the mkcflow model as the 
global deposition rate allowed for an isobaric cooling flow 
across the entire 0.3-3.0 ke'V temperature range. The red¬ 
shift is tied to the best-fit value found in the hottest compo¬ 
nent. This choice does not introduce any uncertainty given 
the strong Ka line complex of the H-like and He-like iron. 

- A set of mkcflow models whose lowest and highest temper¬ 
atures are fixed to cover the 0.3-3.0 keV range with contigu¬ 
ous and not overlapping intervals. Here the normalization 
of each mkcflow component refers to the mass deposition 
rate in the corresponding temperature interval. The upper 
and lower temperatures are frozen to the following values: 
0.3-0.45, 0.45-0.9, 0.9-1.8, and 1.8-3.0 keV. As in the pre¬ 
vious case, a single-temperature mekal model accounts for 
any gas component hotter than 3 keV. The redshift is tied to 
the best-fit found in the hottest component here as well. 

As a fitting method, we used the C-statistics, although we 
also ran our fits with the ;if^-statistics. In the latter case, spec¬ 
tra were grouped with at least a signal-to-noise ratio (S/N) of 4 
in each bin. When running C-statistics, we used unbinned spec¬ 
tra (at least one photon per bin). As a default, we considered 
the energy range 0.5-10 keV both for the MOS and the pn data. 
We quote the best-fit values with Icr error bars on the measured 
value of M, or the Icr upper limit. The same analysis was also 
applied to the Chandra data for a direct comparison. 

4.2. Spectral results on the core region from EPIC MOS and 
pn data 

Before fitting the XMM data, we performed a series of basic 
tests on our spectra. First, we separately fit pn and MOS spectra 
from each exposure using a single mekal model plus the fixed 
AGN component and Galactic absorption. The best-fit tempera¬ 
ture values were all consistent within 1 a, showing that all the 
spectra are broadly consistent with each other. When we added 
a mkcflow component, we noted that the best-fit values of the 
mass deposition rates from pn data were lower than those from 
MOS data in all the different exposures. This shows that there 
may be significant differences between pn and MOS calibration 
that affect the measurements of cold gas components. Therefore, 
we proceeded with our standard analysis strategy, which com¬ 
bines different exposures of the same detector, but fitted pn and 
MOS spectra separately. 

Following our analysis strategy, first we separately fitted the 
MOS and pn data with a single mkcflow model, coupled with a 
mekal model to account for the hot ICM along the line of sight. 
We used both our new XMM-Newton data and the archival data. 
We considered the energy range 0.5-10 keV both for MOS and 
pn data. The XSPEC version used in this work is vl2.8.1. 

First we used a single mkcflow model with !„„„ = 0.3 
keV and Tf„ax = 3.0 ke'V to fit the combined MOS data with 
Cash statistics. The redshift of the cold gas component is linked 
to the best-fit redshift found for the hot component (described 
by a mekal model) zx - 0.593 + 0.002, consistent with the 
optical value (iMcDonald et al.ll2012h . We found that the metal 
abundance of the mkcflow component is not constrained by the 
present data. Therefore we chose to explore the best-fit mass 
deposition rate by varying the cold gas metallicity Zcoid in a 
wide range of values. We set the lowest value for Zcoid equal 
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Fig. 3. Left panel. MOS spectra from the three Obsids with best-fit model obtained for a single mkcflow in the 0.3-3.0 keV 
temperature range. The lower panel shows the residual with respect to the best-fit model. Right panel: components of the best-fit 
model of the MOS spectra (three for each component, shown with dotted lines) obtained for a single mkcflow in the 0.3-3.0 keV 
temperature range. The mkcflow component is the lower thermal component, while the AGN contribution is virtually negligible 
below 2 keV. 
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Fig. 4. Left panel, pn spectra from the three Obsids with best-fit model obtained for a single mkcflow in the 0.3-3.0 keV temperature 
range. The lower panel shows the residual with respect to the best-fit model. Right panel: components of the best-fit model of the 
pn spectra (three for each component, shown with dotted lines) obtained for a single mkcflow in the 0.3-3.0 keV temperature range. 
The mkcflow component is the lower thermal component, while the AGN contribution is virtually negligible below 2 keV. 


to the metallicity measured for the hotter mekal component 
Zhnt - 0.47~^!?ntZf7i (with solar metallicity values measured by 
lAsplund et al.1l2005h . and we set the upper limit to the 2cr upper 
value ~ 1.4Zo found in the inner 80 kpc in the Chandra anal¬ 
ysis. We quote the best-fit values obtained for = 1-0 Zq, 
and add the uncertainty associated to the range of Zcoid as a 
systematic error. We found a global mass deposition rate M - 
620 (-190 -H200)ifa( (-50 -h l50)systMQ yr*'- The intrinsic ab¬ 
sorption of the AGN is found to be NH - (42.0 + 1.5) x 10^“ 
cm“^, in excellent agreement with the value found with the 
Chandra analysis. The best-fit temperature of the hot compo¬ 
nent is kTfiot - (6.5 + 0.3) keV, which also agrees well with 
the emission-weighted temperature in the inner 80 kpc obtained 
from the analysis of the Chandra data. We repeated the fit with 
;^;'^-statistics with spectra binned to at least an S/N= 4 per bin, 
and we found almost identical results, with a reduced = 1.04 
for 663 dof. The binned spectra of the three Obsids with the best- 
fit models are shown in Fig. [2 left panel. We show the different 
components of the best-fit model in the right panel of this figure. 
We note that the contribution of the cold gas in the temperature 
range 0.3-3 keV is about 5% in the 0.5-1 keV energy band, while 


the strongly absorbed AGN has virtually no emission below 2 
keV, also when assuming the maximum uncertainty in the intrin¬ 
sic absorption value. We conclude that the MOS data provide 
a positive detection of an average mass-deposition rate in the 
temperature range 0.3-3.0 keV at the 3cr c.l., with M > \90Mq 
yr“' at the 2cr c.l. after considering the systematic effect asso¬ 
ciated with the unknown metallicity of the cold gas. For com¬ 
pleteness, we repeated our analysis without merging MOS 1 and 
MOS2 spectra, which implies a combined fit of six independent 
spectra. We found similar results and error bars for all the param¬ 
eters, with the best-fit values for the mass deposition rate lower 
by 6.5%. 

Then we ran the multiple mkcflow model on the com¬ 
bined MOS spectra. The constraints on the amount of cold 
gas vary significantly as a function of the temperature range. 
In this case, the metal abundance of the cold gas is linked to 
the value found in the temperature range 1.8-3.0 keV, which is 
Zi. 8 - 3 fev = 0.36^q\^^Zq. The best-fit temperature of the mekal 
component is now kThot - 7.4)^® ^ keV with a metal abundance 
of Zhot - 0.54^g g^Zo. The results for M are shown in Fig. |5] 
left panel. We note that the lowest upper limits (at Icr) are mea- 
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Fig. 5. Left panel: mass deposition rate M as a function of the temperature from the analysis of the MOS data of the Phoenix with 
a multi-component mkcf low model. Error bars are at 1 cr confidence level and arrows represent upper limits at 1 cr (the systematic 
uncertainty is not included here). The red horizontal line shows the best-fit value from the single mkcflow model in the 0.3-3.0 keV 
temperature range, while the shaded area shows the 1 cr uncertainty. Right panel: Same as in the left panel, but for the pn data. 


sured in the temperature bins 0.45 - 0.9 and 0.9 - 1.8 keV. These 
upper limits provide the strongest constraint for the fit with the 
single mkcflow model (see shaded area in Fig.|5|. We also note 
that a clear detection of a mass deposition rate is obtained for 
the 1.8-3.0 keV temperature range. This means that our anal¬ 
ysis of the MOS data with a temperature-dependent mkcflow 
model suggests that the cold gas around and below 1 keV does 
not cool as rapidly as the gas between 2 and 3 keV. This is con¬ 
sistent with the physical conditions often encountered in cool 
cores, where the gas cools down to a temperature floor below 
which it is hard to probe the presence of cooling gas. However, 
in the case of the Phoenix, the gas is observed to cool down 
to a temperature much lower than the ambient temperature TviV, 
as opposed to the classic cool cores where the lowest temper¬ 
ature is ~ 1 /37’vi>- We only obtained upper limits on M for gas 
with temperatures below 2 keV. From comparing the results from 
the multi-temperature mkcflow model with that from the single¬ 
temperature mkcflow, we conclude that the mass deposition rate 
found in the analysis with a single mkcflow model should not be 
taken as representative of the entire cooling flow, since it is an 
emission-weighted value averaged over a temperature range that 
is too wide. Therefore, a temperature-resolved analysis of the 
cool core structure is mandatory to constrain the presence of a 
cooling flow. Nevertheless, the constraints on the mass deposi¬ 
tion rate from MOS data still allow values of several hundreds 
of Mq yr“', which might agree with the star formation rate ob¬ 
served in the BCG. Finally, the same fit on the separate MOS 1 
and MOS2 spectra provides very similar upper limits. In partic¬ 
ular, the most constraining bin at 0.45 < kT < 0.90 keV is left 
unchanged. 

The fit of the pn data with the single mkcflow model en¬ 
countered some difficulties. The best fit was obtained when the 
temperature of the hot component was kT^ot = 5.9 + 0.1 keV, 
which disagrees with the MOS data and with the Chandra re¬ 
sults. In addition, the reduced = 1.38 for 637 degrees of free¬ 
dom, is significantly larger than the value obtained for the MOS 
data. The lower temperature of the hot component is enough to 


account for most of the soft emission, which in the MOS fit was 
associated with the cold gas. As a result, the soft emission that 
can be associated with the cold gas is significantly lower, and the 
3cr upper limit to the mass deposition rate is as low as ~ 24QMq 
yr“'. 

The significant difference between pn and MOS results 
clearly requires a detailed treatment of any cross-calibration un¬ 
certainty, which is beyond the goal of this paper. A first attempt 
to estimate the systematic uncertainty due to cross-calibration 
problems is discussed in Sect. 4.4. Nevertheless, a reliable as¬ 
sumption which may bring the two analyses into better agree¬ 
ment is to set the temperature of the hot gas to the value found 
with Chandra and XMM-Newton MOS, kTkot ~ 6.8 keV. With 
this assumption, we find M = 210 (-80 + S5)stat (-35 -i- 
60)systMe yr Here, the best-fit value was also obtained af¬ 
ter freezing the abundance of the cold gas to Zc„id =1.0 Z©, and 
the systematic error was obtained by varying the abundance of 
the cold gas in the range 0.45-1.4 Zq. The abundance of the hot 
mekal component is Zhot - 0.48 + 0.04 Zq , which agrees very 
well with the MOS data. Finally, the intrinsic absorption of the 
AGN is also well constrained to be NH - (43.6 + 1.2) x 10^^ 
cm“^. The binned spectra of the three Obsids with the best-fit 
models are shown in Fig. |4] left panel. In Fig. |4] right panel, we 
show the different components of the best-fit model. It is pos¬ 
sible to appreciate the lower mkcflow component with respect 
to the right panel of Fig. [3] This means that the analysis of the 
pn data plus the condition kT/to, - 6.8 keV provides a positive 
detection of a significant mass-deposition rate, although about a 
factor of 3 lower than those found in the MOS analysis. 

We briefly comment on the dependence of M on the temper¬ 
ature of the surrounding hot gas. Clearly, the value of M criti¬ 
cally depends on the amplitude of the continuum and therefore 
on the hot gas temperature. As already mentioned, our strategy 
consisted of anchoring the hot gas temperature to the best-fit 
value driven by the continuum in the high-energy range, while 
the colder gas below 3 keV was properly treated by the multi¬ 
temperature mkcflow model. With current data, any other treat- 
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ment of the hot gas component would be highly speculative, 
therefore we did not explore the possible effect of the temper¬ 
ature distribution above 3 keV in the cluster core in more de¬ 
tail. We argue that the best way to deal with this aspect is to 
use highly spatially resolved data that allow avoiding the sur¬ 
rounding hot gas component as best possible and to focus on the 
innermost core region. 

The analysis with a multi-temperature mkcflow model, 
shown in the right panel of Fig. |5] confirms this inconsistency. 
The behavior of M as a function of the temperature is similar, 
with the most constraining upper limits coming from the temper¬ 
ature bins 0.45 - 0.9 and 0.9 - 1.8 keV and a clear detection of 
a high mass-deposition rate in the 1.8-3.0 keV temperature bin. 
The ratio between M in the bin 1.8-3.0keV and the upper limits 
below 1.8 keV also has the same value as in the MOS analysis. 
We therefore qualitatively reach the same conclusion as was ob¬ 
tained from the MOS data analysis with the multi-temperature 
mekal model, but the values M are scaled down by a factor ~ 3. 
Incidentally, if the energy range used for spectral fitting is re¬ 
duced to the 0.7-10.0 keV range, best-fit values and upper limits 
from the pn spectra are found to agree with MOS results. Clearly, 
this simply depends on the fact that the sensitivity to the cold gas 
is significantly reduced by excluding the energy bins below 0.7 
keV, allowing much higher upper limits to the contribution from 
the cold gas. 

To summarize, the analysis of XMM-Newton CCD spectra 
confirms the detection of gas cooling at a rate > IOOOMq yr * in 
the temperature range 1.8-3.0 keV, while it only provides upper 
limits for gas at temperatures kT <1.8 keV. The upper limits to 
the global mass deposition rate below 1.8 keV appear to be con- 
sistent with an SFR of ~ 8 OOM 0 yr ' within Icr, as measured by 
[McDonald et al.l (1201 3h in the BCG, from the MOS data anal¬ 
ysis. On the other hand, the upper limits on M measured from 
pn data at temperatures below 1.8 keV are significantly lower 
(more than 3cr) than the SFR. Therefore, no final conclusion can 
be drawn on the correspondence between the cooling flow and 
the observed star formation rate in the core of the Phoenix clus¬ 
ter. 

4.3. Spectral results from RGS data 

The RGS spectrum was extracted from a strip of roughly 50 arc- 
sec long across the center of the object, a width corresponding to 
90% PSF. Because the RGS is a slitless instrument and the zero- 
point of the wavelength calibration is dependent on the position 
of the detector in the held of view, we used as the source center 
the coordinates RA - 23:44:43.9, DEC - -42:43:12.64, which 
correspond to the position of the central AGN (see Sect. 3). In 
practice, with our choice for the extraction region, we collect 
photons from a box region with dimen sions 50”x 12' cente red on 
the source (see, e.g.. Figs. 2 and 3 in IWerner et al.ll200^ . This 
is the standard choice to maximize the signal from the cluster. 
This choice is motivated also by the fact that the Phoenix cluster 
has a brightness distribution that is similar to a point source from 
the XMM point of view, so that the line widening that is due to 
the spatial extension of the source does not severely affect the 
data. This means that the effective area that is sampled by our 
RGS spectrum is a circle with a radius of 25”, to be compared 
with the circle of 13” used to extract the EPIC spectra. We also 
extracted an RGS spectrum from a narrower region with dimen¬ 
sions 15” X 12', corresponding to 70% of the PSF as opposed to 
the 98% of the PSF achieved with a width of 50”. 

We fitted the first-order spectra between 7 and 27 A, since 
this is the range where the source is higher than the background. 


We usedX^PEC version 12.8.1 and C-statistics on the unbinned 
spectrum. We first applied a single-temperature mekal model 
modified by the Galactic absorption, while the redshift was fixed 
at the optical value. The AGN was modeled as in the fits of the 
EPIC data, but because of the strongly absorbed spectrum, the 
effect of the AGN emission is completely negligible in the RGS 
spectra given the limited wavelength range. Prom the RGS spec¬ 
trum extracted from the larger regions (with a width of 50”), 
we found that a single-temperature model provides a reason¬ 
able fit to the data, with a Cstat/dof — 1978/1952 with a best-fit 
temperature of kT - 6.2 + 0.4 keV with a metal abundance of 
Z = 0.55 + 0.09 Zq. We show in the left panel of Pig.l^the first- 
order combined spectrum with the best-fitting thermal model. In 
the same figure we show the positions of relevant emission lines 
that are clearly identifiable in the spectrum. The set of lines is 
the same seen in the RGS spectrum of the cluster Abell 1835, as 
shown in the comparison plot in the bottom panel of Pig.0 Por 
Abell 1835 we used the same data reduction as performed on 
the Phoenix data on the three observations not affected by flares 
(Obsid 0098010101, 0551830101, and 0551830201). The spec¬ 
trum obtained with our reduction agr ees very well with the spec¬ 
trum shown by ISanders et al.l (l2010h . Based on the RGS spec¬ 
trum of the Phoenix cluster we therefore reach the same conclu¬ 
sions as were obtained for A1835: no line emission from ioniza¬ 
tion states below Pe XXIII is seen above 12A, and no evidence 
for gas cooling below ~ 3 keV is found. 

We then proceeded to obtain a measurement of the mass de¬ 
position rate from the RGS data by using the mekal+mkcflow 
model with abundances constrained to be the same in both com¬ 
ponents. As for the EPIC data analysis, the temperature range of 
the mkcflow component was fixed to 0.3 - 3.0 keV. We found no 
statistically significant improvement with respect to the single¬ 
temperature model (Citof/dof = 1977/1951). The best-fit ambient 
hot temperature was kThot — (6.3 + 0.5) keV with a metal abun¬ 
dance of Z = 0.43 + 0.09 Zq. The best-fit value for the mass 
deposition rate was M - 122/j22 yr *■ The 90% upper limit 
on the mass deposition rate was 682 M© yr^'. 

More complicated models, such as the one including a set of 
mkcflow models, are not constrained by the RGS data. However, 
we were able to obtain constraints on the abundance of metals 
other than iron using a vmekal model. We found Z^g = (0.53 + 
0.24)Zo,Zs, = (0.66 + 0.15) Zo, andZj = (0.33 + 0.25) Zq. The 
use of vmekal provides a slightly lower upper limit on the mass 
deposition rate: M - 103/jq3 Mq yr“*, with a 90% upper limit 
of 590 Mq yr '. Best-fit temperatures are left unchanged. 

We applied the same analysis to the RGS spectrum extracted 
from a smaller region with a width of 15”. We obtained very 
similar results, with small differences, the largest being a lower 
value of about 1 cr of the ambient hot temperature, which was 
meaured to be kThot - 5.7 /q3 keV. As for the mass deposition 
rate, we found a 1 cr upper limit of M = 470 M© yr“', with a 90% 
upper limit at 744 M© yr"'. Therefore, the RGS spectral analysis 
provides similar results for the two different extraction regions. 

To summarize, the analysis of the RGS data only provided 
upper limits to the global mass deposition rate, in agreement 
with the values found with the analysis of the EPIC pn data, 
and also consistent with those found with EPIC MOS within 
less than 2cr . The better agreement with EPIC pn analysis is 
expected on the basis of the cross-calibration status of XMM in¬ 
struments: the agreement between EPIC pn and RGS is currently 
within a few percent as a result of the combined effect of the 
RGS c ontamination mo del plus the improved pn redistribution 
model (IStuhlingeill20ir)h . 
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At the same time, we remark that the comparison between 
the values of M obtained from the two RGS extraction regions 
and the EPIC pn extraction region allows us to conclude that our 
choice of the extraction radius for EPIC data was adequate for 
sampling the possible cold gas emission. The simple fact that M 
from RGS sampled in the region twice larger than that used for 
EPIC data is not significantly larger than the value found in the 
pn data confirms that no significant emission from cold gas can 
be found outside the 13” radius. 


4.4. Discussion of the systematics in XMM-Newton data 
anaiysis 

The discrepancy among the results obtained with MOS and pn 
can be ascribed to the uncertain cross-calibration of the effective 
area among the XMM-Newton CCD, while there are no other 
sources of statistical noise that can account for a significant frac¬ 
tion of this discrepancy. This difference is persistent despite the 
recent model for the MOS contamination intr oduced in the SAS 
release vl3.5.0 (see lSembav & Saxtonll201^ . Since the contri¬ 
bution of the cold gas below 3 keV to the total emission in the 
extracted region is about 3% in the 0.5-2.0 keV energy range and 
about 5% in the 0.5-1.0 keV range, an uncertainty in the effec¬ 
tive area of the same order may severely affect the measurement 
of the cold gas. 

A first way to estimate the calibration uncertainty is to com¬ 
pare our results with the results obtained by applying the cor¬ 
rection CORRAREA to t he response area of the EPIC detec¬ 
tors based on the results of iRead et al.l(l2014l) . The CORRAREA 
calibration is based on an extensive cross-calibration study 
of 46 non-piled-up sources ext racted from the 2X MM EPIC 
Serendipitous Source Catalogue (IWatson et al.ll200^ . This phe¬ 
nomenological correction is meant to bring into agreement the 
broadband fluxes measured by EPIC-MOS and EPIC-pn. We 
obtained an estimate of the uncertainty on M values associated 
to uncertainties in the EPIC response areas by comparing the 
best-fit values obtained with and without the CORRAREA cor¬ 
rection. The best-fit values of M are 1% and 5% higher when 
applying the CORRAREA correction for MOS and pn data, re¬ 
spectively. This is far from the factor of 3 needed to bring M 
from MOS and pn into agreement. Therefore we conclude that 
the measurement of cold gas in the spectra of the Phoenix cluster 
is significantly affected by the c alibration of the E PIC detectors 
at a level beyond that probed bv iRead et al.l (l2014h . 

We also included the energy range 0.3-0.5 keV to investigate 
possible effects associated with this low-energy band. Although 
the calibration is even more uncertain in this range, the inclusion 
of the lowest energies may help to increase the relative contribu¬ 
tion of the cold gas. However, when we repeated our fits on the 
0.3-10 keV energy range, we did not find any significant differ¬ 
ence with respect to the results described in Sect. 4.2. 

We also explored the possibility that data taken in different 
epochs may have a different calibration. However, when we ex¬ 
cluded the shortest Obsid, which was acquired two years earlier 
than most of our data on the Phoenix, the best-fit values were 
only affected by a negligible amount. 

A main source of uncertainty in the soft band is the Galactic 
column density. The possible presence of unnoticed fluctua¬ 
tions in the Galactic neutral hydr ogen column densitie s on 
scales smaller than the resolution of iKalberla et al. ( 20051) ma y 
reach 10-40% on scales of ~ 1 arcmin (iBarnes & Nulsenll2003h . 
Therefore, we considered a systematic uncertainty in our best-fit 
values of M assuming a maximum variation of NHcai by 40%. 
The typical effect is that the best fit is obtained for the high- 



Fig.7. Mass deposition rate M as a function of the temperature 
from the analysis of the ACIS-I Chandra data of Phoenix within 
30 kpc from the center. The AGN emission has been removed 
(inner 1.5 arcsec). The results have been obtained with a multi¬ 
temperature mkcflow plus a mekal model for the hot compo¬ 
nent. Arrows refer to 1 cr upper limits. The shaded area shows 
the 1 cr upper limit based on the analysis of a single-temperature 
mkcflow model in the 0.3-3.0 keV temperature range plus a 
mekal model for the hot component. 


est allowed value of NHcai, and this implies that much more 
cold gas can be allocated. Specifically, if we allow the Galactic 
absorption to vary up to NHcai - 2.13 x 10^° cm“^, we find 
M - 900(-390 -I- 110)j,a,(-130 h- 250)jj,j,Mo yr“* from the 
analysis of the MOS data with a single mkcflow model in the 
temperature range 0.3-3.0 keV. This v alue is 45% higher th an 
that obtained with the NHcai value of IKalberla et al.l (l2005h at 
the position of the cluster. A similar effect is found for the other 
fits, including those with the multi-component mkcflow model. 
If, on the other hand, we leave NHcai free to vary, we find best-fit 
values of NHbestfit - (2.2 + 0.8) x 10^’’ cm“^, which is very close 
to the upper bound we assumed for NHcai- The analysis of the 
EPIC pn data, with the same values of NHcai, provides M values 
2.3 and 3 times higher for NHcai = 2.13 and 2.2 x 10^“ cm“^, re¬ 
spectively, which considerably reduces the discrepancy between 
the MOS and pn analysis. However, when the Galactic absorp¬ 
tion is left free, the best-fit value is NHcai = (0.8 + 0.2) x 10^° 
cm“^, which in turn provides much lower values of M than the 
standard analysis. Therefore, we conclude that there are no hints 
for a plausible variation of NHcai that can significantly change 
our results, including the discrepancy between the best-fit values 
of M found between MOS and pn data. 

4.5. Comparison with spectrai resuits on the core region 
from Chandra data 

We repeated the same fits on the Chandra data available in the 
archive as of December 2014. The photometry in the inner 30 
kpc, corresponding to 4.5 arcsec, excluding the central AGN, 
amounts to 1510 net counts in the 0.5 - 7.0 keV energy band. 
This is usually not sufficient to provide robust constraints on the 
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Phoenix RGS data 



A1835 end Phoenix data 



Fig. 6. Top panel, combined RGS 1 and 2 spectrum for the Phoenix cluster and best-fit single-temperature model. The data have 
been rebinned to have a signal-to-noise ratio of 10 and divided by the effective area of the instrument only for displaying purpose. 
Relevant lines are also labeled. Bottom panel. The Phoenix spectrum (in red) is compared with the RGS spectrum of Abell 1835 (in 
black). We plot the fluxed spectra obtained with RGSFLUXER. Similar results are obtained by plotting the unfolded spectrum in 
XSPEC (the green line in the Phoenix spectrum) wit the only difference set by the binning condition S IN - 10. The best-fit thermal 
model of the Phoenix is also plotted as a blue line. We also plot relevant lines clearly seen in the spectra: the same sets of lines are 
present in both spectra in the overlapping wavelength range. 


amount of cold gas in nearby clusters, and it would be even less 
effective at the redshift of Phoenix. Nevertheless, we performed 
our spectral analysis as for the XMM-Newton data. The single¬ 
temperature mkeflow model in the temperature range 0.3-3.0 
keV provides a 1 cr upper limit of 630 (-85 -H 60) sysi Mq yr“* or 
a 2cr upper limit of 1130 (-150 + 1100)svs, Mq yr“*, where the 
systematic uncertainty corresponds to 0.45 < Zcoid < I.TZq. The 
hot gas temperature is kTi,oi - keV, within Icr from the 

value found in XMM-Newton analysis and in the overall tem¬ 
perature profile of the Chandra data. We recall that we did not 


fit the AGN since its emission was removed from the spectrum 
thanks to the Chandra angular resolution. 

Despite the low S/N of the Chandra data (almost two or¬ 
ders of magnitude fewer photons than for the combined XMM- 
Newton spectra), and the high redshift of the Phoenix cluster, we 
also performed the fit with a multi-component mkeflow model. 
Temperature and metallicity of the hot gas component were the 
same as in the previous fit, while Zcoid = 1.0 Zq. As expected, 
we obtained little additional information. The results are sum¬ 
marized in Fig.|7] where we can conclude that the strongest con- 
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Fig. 8. Spectral energy distribution of the BCG of Phoenix in the 
FIR from the PACS (100 and 160 ^lm) and SPIRE (250, 350,500 
^m) instruments onboard the Herschel Space Observatory. The 
red dotted line and the blue dashed line show the best fit to the 
SED for the AGN and starburst template, respectively. 


straints on M mostly comes from the gas between 0.9 and 1.8 
keV. This means that Chandra upper limits are consistent with 
a mass deposition rate of about ~ IOOOMq yr“'. We remark, 
however, that deeper Chandra data are expected to provide much 
stronger constraints. 

5. Exploring the connection between SFR and 
global M 

We inspected available EIR data on the Phoenix cluster to bet¬ 
ter constrain the SER. The BCG has been observed in the 
EIR with the PACS (100 and 160 /urn) and SPIRE (250, 350, 
500 um) instrument s onboard the Herschel Space Observatory 
(iPilbratt et al.ll2010h . Herschel data bracket the critical peak of 
EIR emission of high-redshift galaxies, providing a direct, unbi¬ 
ased measurement of the dust-obscured star-formation. We per¬ 
formed aperture photometry assuming an aperture radius of 6 
and 9 arcsec for 100 and 160 /um, respectively, while for SPIRE 
we used PSF fitting. The five data points from Herschel instru¬ 
ments provide the SED shown in Eig. |8] 

We investigated the contribution of the AGN compo- 
nent to the FIR em ission using the program DECOMPIR 
dMullanev et al.ll201 ill . an SED model-fitting software that aims 
to separate the AGN from the host star-forming (SE) galaxy. The 
AGN component is an empirical model based on observations of 
local AGNs, whereas the five starburst models were developed 
to represent a typical range of SED types, with an extrapolation 
beyond 100 /im using a gray body with emissivity jd fixed to 1.5. 
The best-fitting model obtained with DECOMPIR confirms at 
least a 50% contribution of the AGN to the total SED flux, as 
shown in Eig. [8] The AGN clearly dominates the EIR emission. 
The IR luminosity fo r the starburst compon ent (best described 
by the model SB5 in iMullanev et al.ll201 ih is measured to be 
L/r = 3.1 X 10^^ Lq. Assuming a Salpeter IME, this luminos¬ 
ity corresponds to a star formation rate of 530 Mq yr“' with a 
typical error of 15%, which is lower than the value found by 
iMcDonald et al.l (120131) at 3 cr c.l. The total luminosity given by 
the best fit is Ljr - 2.09 x IO'^Lq, with the AGN component 
contributing by 86%. Taken at face value, this revised estimate 
of the SFR is in very good agreement with the mass deposition 
rate found with the analysis of the EPIC MOS data, while is still 


inconsistent at about 3cr c.l. with the mass deposition rate found 
with EPIC pn analysis. 

6. Conclusions 

We analyzed the X-ray data taken with a 220 ks exposure of 
XMM-Newton on the Phoenix cluster. We focused on cold gas 
in the core by selecting a circle of 13.5 arcsec centered on the 
BCG in the XMM image and a strip with a width correspond¬ 
ing to 90% of the PSP in the RGS data. Our immediate goal 
was to constrain the actual mass-deposition rate associated with 
the residual cooling flow, with the aim of understanding whether 
this may be directly linked to the current massive starburst ob¬ 
served in the BCG. We combined XMM-Newton data with shal¬ 
low Chandra data, particularly to model the hard emission of 
the central absorbed AGN, which cannot be removed from the 
XMM-Newton data alone. Our results are summarized as fol¬ 
lows: 

- We measured an average mass-deposition rate of M - 
620 (-190 -H 200)„a, (-50 h- 150),^,, Mq yr^' and M = 
210 (-80 +85W (-35 +60V,MG,yr-'inthe0.3-3.0keV 
temperature range from the analysis of the MOS and pn data, 
respectively. These values are dominated by the cold gas in 
the energy range 1.8-3.0 keV, while only upper limits can be 
obtained at temperatures below 1.8 keV. The upper limits to 
the global mass deposition rate below 1.8 keV appear to be 
consiste nt with an SFR of ~ SOO Mp yr ', as measured in the 
BCG bv lMcDonaldet al.l (l2013l) . from the EPIC MOS data 
analysis within Icr, while the upper limits on M measured 
from EPIC pn data at temperatures below 1.8 keV are sig- 
nificantly lower (more than 3cr) than the SFR measured in 
iMcDonaldetaP (l20Tl . 

- Considering the temperature range 0.3-1.8 keV, we found 
that MOS data analysis is consistent with M ~ 1000 within 
Icr, while the pn data proyide M < 400Mq yr“' at 3cr c.l. 

- Since the discrepancy between the MOS and pn data anal¬ 
yses cannot be explained on the basis of currently known 
cross-calibration uncertainties between the two instruments 
nor of other sources of statistical noise, we argue that addi¬ 
tional calibration problems between EPIC instruments still 
need to be understood and properly treated. The contribu¬ 
tion of the cold gas with an ayerage mass-deposition rate of 
~ 600Mg yr ' in the EPIC data is about 5% in the energy 
range 0.5-1.0 keV for our extraction region. This implies that 
any calibration uncertainty on the same order strongly affects 
the data. This conclusion is yalid in the framework of the iso- 
baric cooling model that we assumed here. At present, we are 
unable to discuss whether the assumptions of different phys¬ 
ical models for describing the cooling of the gas can change 
the picture and mitigate the difference between MOS and pn 
analysis. 

- No line emission from ionization states below Ee XXIII is 
seen aboye 12A in the RGS spectrum, and the amount of gas 
cooling below ~ 3 keV has a formal best-fit yalue for the 
mass deposition rate of M - 122^j22 Mq yr“'. This result 
was confirmed by a direct comparison of the Phoenix RGS 
spectrum with A1835 in the oyerlapping spectral range. This 
means that the mass deposition rate f r om R GS analysis is 
lower than the SFR of lMcDonald et al] (120131) in the BCG at 
the 2cr c.l. 

- Current Chandra data (from a short exposure of ~ 10 ks) 
agree with our XMM-Newton analysis, but do not proyide 
meaningful constraints on M. Deeper Chandra data, thanks 


12 


























Tozzi et al.: XMM-Newton observation of the Phoenix cluster 


to the high angular resolution, are expected to provide tighter 
constraints to the global mass deposition rate. 

- A careful analysis of the FIR SED based on Herschel data 
provided a value for the SFR in the BCG of 530 Mq yr“* 
with an uncertainty of 15%. This revised estimate of the SFR 
agrees very well with the M from EPIC MOS data and is 
consistent within Icr with M from the RGS analysis, while 
still inconsistent at more than 3cr with the mass deposition 
rate found in EPIC pn data. Our revised SFR therefore does 
not significantly change the comparison between SFR and 
the global M in the Phoenix cluster. 

To summarize, the range of M allowed by XMM-Newton 
data from EPIC MOS is consistent with the SFR observed in the 
BCG, while EPIC pn and RGS data analyses suggest a mass de¬ 
position rate of about ~ SFR/3, and the derived upper limit is 
inconsistent with the observed SFR at least at the 2cr level for 
RGS, and more than 3cr for EPIC-pn. As a consequence, our 
results do not provide a final answer on the possible agreement 
between the mass deposition rate of isobaric cooling gas in the 
core and the observed SFR in the BCG. As recently shown by 
Molendi et al. (in preparation), M is often measured to be sig¬ 
nificantly lower than the global SFR in central cluster regions 
for several strong cool-core clusters. These findings suggest that 
cooling flows may be short-lived episodes, efficient in building 
the cold mass reservoir that, on a different timescale and possi¬ 
bly with some delay, triggers the star formation rate in the BCG. 
To investigate whether this also occurs in the Phoenix cluster, or 
whether the Phoenix cluster actually hosts the highest cooling 
flow observed so far with M ^ S FR, we must wait for a deep, 
high-resolution, spatially resolved X-ray spectral analysis to re¬ 
move the stronger emission from the surrounding hot gas as best 
possible. 
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